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We introduce coupled channels of Bethe-Salpeter mesons 
both in the boundstate equation for mesons and in the mass 
gap equation for chiral symmetry. Consistency is insured 
by the Ward Identities for axial currents, which preserve the 
Goldstone boson nature of the pion and prevents a systematic 
shift of the hadron spectrum. We study the decay of a scalar 
meson coupled to a pair of pseudoscalars. We also show that 
coupled channels reduce the breaking of chiral symmetry, with 
the same Feynman diagrams that appear in the coupling of a 
scalar meson to a pair of pseudoscalar mesons. Exact calcu- 
lations are performed in a particular confining quark model, 
where we find that the groundstate I — 0, 3 Po qq meson is 
the /o(980) with a partial decay width of 40MeV. We also 
find a 30% reduction of the chiral condensate due to coupled 
channels. 
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I. INTRODUCTION 



The Scalar puzzle 



The scalar mesons form perhaps the most puzzling 
family in hadronic physics. The first puzzling fact con- 
cerns the experimental errors in the partial decay widths, 
the decay widths and even the masses. The lightest 
scalar, /o(400 — 1200), has a poorly determined mass. 
The confidence on the decay widths of the /o(980) and 
a (980) has also decreased Q strongly since 1994. 

The other puzzling argument concerns the matching 
of the the nine observed states with simple SU (3) / qq 
states. One would expect four different towers of states 
corresponding to the two / = /o (which are not de- 
generate for instance because the quark s has a clearly 
larger mass than the quarks u, d) and to the I = 1 ao , 
and the 1 — 1/2 Kq. A short glance at Fig. [l]is sufficient 
to discard the single light and extremely broad state, the 
Jo (400 — 1200) as a simple member of this family. Then 
for the groundstates we could ascribe the narrowest states 
which are respectively the a (980), i^(1430), / (980) 
and /o(1500), and for the radial excited states we could 
respectively ascribe the a (1450), #£(1950), / (1370) 
and Jo (2200). However the decay widths T of the ground- 
state scalars are narrower than expected when compared 
with other resonances decaying in the same pseudoscalar 
pairs but with higher angular momentum, except for the 
only precisely measured one, the ATq(1430). Moreover 



the breaking of SU(3)f due to the m u ~ md « m s 
mass difference is nearly the double than expected when 
compared with the vector meson family and with most 
baryons. This is only comparable with the splittings in 
the pseudoscalar family. 




FIG. 1. We represent, with rectangles of height T, the 
experimentally observed scalar resonances according to the 
Review of Particle Properties, and mark the unconfirmed ones 
with a " ?" . In a naive quark model, the ao and Kq channels 
should show half of the resonances of the /o channel. 

The scalar family is also the most interesting place to 
search for the lightest (S-wave) non qq states. There are 
several theoretical candidates to extra states which may 
be found. The lightest glueball, which is expected from 
QCD, should be a scalar ||. The model of Isgur and 
Weinstein et al (|] suggests that the narrowest scalars 
are meson-meson molecules. The One Meson Exchange 
Potential models for the NN interaction usually postu- 
lates a scalar meson a with a light M ~ 0.5GeV. In 
strongly coupled effective meson models ||,^] , extra poles 
appear in the S matrix when couplings are large. These 
meson models turn out to be the most successful models 
ones so far, explaining with complex nonlinear effects not 
only the narrow ao(980) and /o(980) which are due to the 
vicinity of the KK threshold, but also the very wide and 
light / (400- 1200). 



The relevance of chiral symmetry breaking 

Chiral symmetry breaking is important for the study 
of scalar mesons and their decays for several reasons. Un- 
like the vector, axial and tensor mesons, the scalar and 
pseudoscalar mesons are mixed by the chiral rotations, 



ipip — > co8(6)tpip + i sin(9)ip"f5ip 



(1) 



1 



thus scalars and pseudoscalars are particularly sensitive 
to the chiral symmetry. The very small mass of the 
pseudoscalars is usually explained with chiral symme- 
try breaking. Moreover, since scalars decay essentially 
in pseudoscalars, the pseudoscalar mass is important for 
the scalar decay. Thus we expect that not only the pseu- 
doscalar mesons but also the scalars must contain the 
signature of the breaking of chiral symmetry. Inversely, 
the breaking of chiral symmetry is generated from the 
trivial vacuum by scalar condensation, and we also ex- 
pect that the scalar properties should affect the breaking 
of chiral symmetry. 

These effect have been studied in phenomenological 
meson sigma models, see for instance (|| . At the more mi- 
croscopic level of quarks, Dynamical Spontaneous Chiral 
Symmetry Breaking (x^B) has been worked out in the 
past with several different quark-quark effective interac- 
tions, at the same level as Bardeen Cooper and Schrief- 
fer (BCS) did || for Superconductivity. Sinc e Nambu 
and Jona-Lasinio (NJL) and until recently |^|,^Jthe mass 
gap equation in chiral physics has been so far of the BCS 
type, j?| including only the first order contribution from 
the quark-quark interaction. In this case the quark con- 
densate consists of scalar 3 P quark antiquark pairs. 
At the BCS level, which is very consistent, it is possible 
to derive the mass gap equation in several different but 
equivalent methods. The pseudoscalar meson properties 
have been studied in great detail and the full meson spec- 
trum has been also been calculated in the literature . 
However the BCS approach is not exact in the case where 
coupled channels of mesons are included. 



Coupled Channels of Mesonic qq pairs 

In the case of weak coupled channel effects, it would 
be acceptable to start from boundstates obtained at the 
BCS level and couple them with the help of the anni- 
hilation diagrams of Fig. ||, without changing the mass 
gap equation. In this sense we started some years ago 
to develop a program jl2| to study the coupled channel 
effects in quark models with chiral symmetry breaking. 




(a) (b) (c) 

FIG. 2. We show the quartic diagrams which may con- 
tribute to the boundstate equation. In the strong coupling 
BCS, (a) is included in the self energy and (b) is included 
in the interaction kernel. Diagram (c) which creates or an- 
nihilates quark-antiquark pairs is only used in what we call 
beyond BCS. 

The first result of our program was to reproduce jl^ ] 
at the BCS level the strong decay of the vector meson 
p (and of the < $). This has been studied by other au- 
thors recently |L3fl. Later we extended our program to 



the nucleon interactions and had good results |12| in the 
K n s-wave scattering, the F^^n and Faoa derivative 
couplings and the N N short range interaction. 

However there is a recent trend in the literature to 
reevaluate coupled channel effects in many hadronic phe- 
nomena. Some years ago they were not supposed to ac- 
count for more than 10% of a hadron mass but presently 
they are supposed to contribute with a negative mass 
shift of the order of 50% of the bare mass Jl^Irl . 

Moreover it is possible to prove that the vacuum so- 
lution of the BCS mass gap equation is not exact when 
coupled channels are included. Suppose that the mass 
gap equation for chiral symmetry breaking was solved at 
the BCS level, i.e. without including the coupled chan- 
nels, then a bare pion with vanishing bare mass would 
be found. If the coupled channels were then included, 
at posteriori, in the bound state equation then the pion 
mass would be the sum of the small bare mass plus a 
mass shift and thus would have a resulting nonvanishing 
mass, which implies that the pion lost its Goldstone bo- 
son nature. This result is unavoidable, and can be proved 
variationally. Thus the mass gap should be solved beyond 
BCS, especially when scalar mesons are studied. 

Solving this problem is a corner stone of our program. 
It amounts to join the BCS mechanism with the mean 
field expansion of effective mesons. The logical path of 
the method we will follow is illustrated in Fig. ||. Self 
consistency is insured by Ward Identities, and the scalar- 
pseudoscalar coupling will turn out to be crucial for this 
development. The reward of solving chiral symmetry 
breaking with coupled channels is a 7r with a vanishing 
positive mass in the chiral limit complying with all the 
theorems of PCAC, and a tower of resonances (including 
the scalar meson resonances) above the it with higher 
masses due to radial, angular, or spin excitations. In 
particular the resonances also have an imaginary compo- 
nent of the mass =lf- that describes the decay width into 
the open channels. 



BCS: 




beyond 
BCS: 

FIG. 3. We illustrate the principle which is followed in this 
paper to include the coupled channels in the mass gap equa- 
tion 



The paper 

The aim of this paper is to study at the quark level 
some meson decays which were studied in the literature 
without including directly the full quark contributions 
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We specialize in the groundstate fo — > ir ir decay. 
We also study the effect of the meson coupled channels 
on the quark condensed vacuum. 

The remaining of this paper is organized as follows. In 
section II we review the scalar masses and decays at the 
BCS level. This includes the choice of an effective inter- 
action for quarks, the mass gap equation and the Bethe 
Salpeter equation at the BCS level, the scalar coupling 
to pseudoscalars, and the scalar decays width. In Sec- 
tion III we produce a finite extension beyond BCS for a 
class of confining effective interactions, derive the mass 
gap equation with coupled channels, link it to the scalar- 
pseudoscalars coupling, and solve the mass gap equation. 
Results are shown in Section IV together with their dis- 
cussion. We also include 4 appendices. 



II. THE BCS LEVEL FOR A PARTICULAR 
FORMALISM 

A. The choice of an effective interaction 

The quantitative results of this paper will be obtained 
with a particular chiral invariant strong potential which 
is an extended version of the N JL potential . For the 
study of dynamical x^B it is crucial to have a closed 
model where calculations can be carried until the end, 
because precise cancellations occur. At this point we 
abandon the explicit SU(3) gauge invariance. We start by 
introducing a class [Tl]Jl7| | of Dirac quark Hamiltonians 
which arc, in the limit of massless quarks, explicitly chiral 
invariant, 



H = I d 3 x 



ip^(x) (mo/3 — ia. ■ V) ip(x) + 



iP(x)Slrt(x) d 4 y Vt(x - y) jp(yWil>{y) 



(2) 



The quark-quark interaction, is an effective color depen- 
dent 2-body interaction. In eq. @) the operators Qi in- 
clude both the color Gell-Mann matrices and the Dirac 
matrices. The sum in the Dirac matrices must be chiral 
invariant. Because the Gell-Mann matrices are traceless 
there will be no tadpoles in this scheme. In Hadronic 
Physics the effective interaction should be simultane- 
ously color confining, in the Minkowsky space, local, and 
Lorentz invariant. However no interaction which com- 
plies with all these constraints has yet been used to study 
chiral symmetry breaking. 

The models used in the literature divide essentially in 
two classes, in particular the more popular one springs 
directly from N JL , is Euclidean and due to the struc- 
ture of the interaction has usually no analytic continua- 
tion to the Minkowsky space and lacks confinement. This 
class has been extended in many different directions. For 
instance finite size boundstates were included with the 



Global Color Model of ref. |l8[ , and a sophisticated in- 
teraction with a general tensor structure, an almost linear 
long range and pertubative short range is found in ref. 

& ___ 

The other class of models |lG,nilE0| has the single 
drawback of using an instantaneous potential (except for 
Lorentz invariant extensions pl[ of this class). But it has 
the advantage of being confining which allows to study 
the whole [|n],^2|jl|,|22]-||| hadron spectrum. This ap- 
proximation also has the advantage to allow a straight- 
forward application to low energy nuclear physics. For 
the theoretical foundations of these models, including 
the connection to both pertubative and nonpertubative 
QCD, see[||. 

Thus we choose to calculate the quantitative results of 
this paper within the second class of Nambu and Jona- 
Lasinio potentials. We use a simple model which is in 
very good agreement with the experiments in what con- 
cerns the hadronic spectroscopy, the decays of the vec- 
tor mesons p and <&, the coupling of a n to a N or A 
and the N N short range interaction moreover it 
supports that chiral symmetry breaking is very stable in 
the presence of Nuclear Matter . While confinement 
is an essential physical aspect of the model, the instan- 
taneous approximation simplifies drastically the energy 
dependence of the interaction, and allows to work in a 
framework which is familiar to Schrodinger' s equations. 
The 2-body potential for Dirac quarks is, 



-3 A A 



. ., 2 [70 <8> 70 (k 3 r 2 -U)+a^- ®jk 3 r 2 ] 6(f) (3) 
and the Fourier transform of the potential is, 

~ —3 A A 

fliVi(k) <g> fi ; = — - ■ ®- [70 <8> 70 (~K 3 A k - U) 

+arf-&?(-K$A k )] (2n) 3 5 3 {k) , (4) 

where we dropped the sum in color and Dirac indices. 
The factor —3/4 simplifies the color contribution for color 
singlets. 



V(k) 




FIG. 4. 



-Ue- 



is an example of a potential which 



tends to —U + k^r 2 in the limit of infinite U. We illustrate 
this in the case where Kq = 1, U = 100. 
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The 70 • 70 term in the potential is the limit of a series 
of attractive potentials, see Fig. |J, and U is an arbitrarily 
large infrared constant. The infinite U reappears in the 
self energy and in color singlet channels this cancels the 
infinitely attractive potential. Any colored state will have 
a mass proportional to U and will thus be confined, see 
appendix A. 

The choice of an harmonic potential is not crucial, a 
linear |Q or funnel potential has also been used, but a 
quadratic form is simpler. In the case of light quarks 
the current quark mass m is almost vanishing and it 
essentially affects the family of the tt which is a quasi 
Goldstone boson. 

The 7 • 7 term is introduced in order to have Lorentz 
invariant pions which are relativistic in the scalar decay. 
Clearly the wrong result of a simple 70 <8> 70 instanta- 
neous interaction is the constant f K which is quite small 
and is not Lorentz invariant. The hope to cure / w with 
covariant extensions of the model turned out to fail since 
they merely plj increased f„ in 30%. The value for the 
parameter a which renders the pion Lorentz invariant is 
a = —.18, and for most calculations (except for f n which 
is increased by 300%, see subsection D) this a yields a 
result comparable to the one we had for a = 0. The small 
a suggests that the cure of f n may be related with the 
pertubative short range pEf quark-quark interaction . 

Once a is fixed, this model has the single scale of the os- 
cillator parameter Kq. The simplest adimensional units 
of Kq = 1 will be used from now on in computations. 
When comparing with experiments we will rescale K to 
the value of Kq — 330 M eV which gives the best over- 
all fit of the meson spectrum. We will also work in the 
momentum representation and drop the tilde ~ from the 
potential. 



B. Chiral Symmetry Breaking at the BCS level with 
quarks and antiquarks 

At the BCS level and with a color dependent inter- 
action the Schwingcr Dyson equation for the quark self 
energy (which is also the mass gap equation) is, 



1 = Sq 1 - 



S- 1 = -i p- E 



(5) 



where the full (up to the approximation which is chosen) 
propagator S is denoted as usual by , . The subindex 
is reserved for the free functions. The effective quark- 
quark 2-body interaction of eq. (Q), a chiral invariant 
and color dependent interaction, is represented with a 

dotted line The Bethe Salpeter equation for the 

vertex is related by the Ward Identities with eq.(^), see 
Appendice C, 



= r n 



d A k 



r(p, q) = r - 

T(p + k,q + k)S(q + k)Q 



■V(fc)fiS(p + AO 



(6) 



where the full vertex T M is denoted by •= . When eq. (g) 
is iterated, we find that it includes, the Bethe-Salpeter 
ladder, which will be represented diagrammatically by a 
box with 4 emerging lines 



=0= 



+ 



+ . 



(7) 



where the ladder represents the mesons, see Appendice 
C. In this way the quark propagator, the vertices and the 
mesons are intertwined. 

In this case of an instantaneous interaction, it is conve- 
nient to substitute the Dirac fermions in terms of Weyl 
fermions, in order to find the hadron spectrum. The 
Dirac propagator can be decomposed in a quark propa- 
gator and an antiquark propagator, moving both forward 
in time. 



SDirac(ko, k) 



jk — m + it 

i 1 + ^P + ^a-k 



ko — E + ie 
i 



P 



l-f/?--§q-fc 

k Q - E + ie 2 



p (8) 



It is convenient to use |TT],|lo| the quark energy projectors, 
A+ = ~ (l + S/3 + d ■ a) = U M 



A" 



\ (l - SP - Ck ■ a) 



(9) 



where S = sin(ip) = ^ , C = cos(ip) = -g and ip is 
a chiral angle which in the non condensed case is equal 
to arctan ^ , (mo is the current mass of the quark) but 
is not determined from the onset when chiral symmetry 
breaking occurs. In this case the physical quark mass 
is a variational function m — m(k) which is determined 
by the mass gap equation. This is equivalent to use the 
chiral angle if = <p(k) as the variational function. In 
Fig. U we show examples of non trivial solutions for the 
function (p(k). 

The energy projectors can be decomposed in the quark 
spinor u(k) and in the antiquark spinor v(k), 



u s (k) = 
«.(k) = 



A" 



1+5 
2 



=«.(0) 



1+S l-S 



-k ■ a 



A- 



1+5 
2 



=».(0) 



1+S l-S 



-k- 



-icr 2 7 5 M*(k) 



u,(0) 

v s (0) 
(10) 
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And finally the Dirac quark propagator is decomposed 
in, 

S Dir ac{w,k) = u(k)S q {po,k)u^(k)/3 

-vHk)S 9 (- Po ,-k)v(k)(3 , (11) 
where the quark and antiquark Weyl propagators are. 



S q (w,k) = Sq(w, k) 



E{k) 



(12) 



The quark and antiquark formalism is convenient to cal- 
culate the hadron spectroscopy. With Weyl propagators 
the BS equation simplifies into the Salpeter equation, in 
a form which is as close as possible to the more intuitive 
Schrodinger equation. In the Feynman rules with Weyl 
propagators, we choose to redefine the vertices of the ef- 
fective potential which now include the spinors , u , 
and v. The " — " sign which affects the antiquark propaga- 
tor in cq. (|lT|) could also be included in the vertices with 
v', but we prefer to recover the equivalent rules which 
are common to nonrelativistic field theory. This "— " 
sign together with the one from the fermion loops will be 
included in the antiquark vertex and in diagrams with 
quark exchange or with antiquark exchange. The Dirac 
vertex 70 is now replaced by u'u , u'v , v^u or — v'v when 
the vertex is respectively connected to a quark, a pair cre- 
ation, a pair annihilation or an antiquark; and the Dirac 
vertex 7 is respectively replaced by u^au ,u^dv , v^au or 
— v'av. We choose the graphical notation for the Weyl 
propagators of quarks and antiquarks, 

w.k 

S D irac{w, k) = D , 



Sq{w, k) 



Sq(—W, -k) 



(13) 



where the Diagrams using the Feynman rules correspond- 
ing to the Dirac fermion propagators will have a subindcx 
£> in the remaining of the paper. In the case of the Weyl 
propagators (which will be used more often the the Dirac 
propagators) the quark will be represented with an arrow 
pointing to the left while the arrow pointing to the right 
represents an antiquark (both move forward in the time 
direction). 



C. The BCS mass gap equation and the quark 
energy 

Here we derive the mass gap equation, and the quark 
dispersion relation, replacing the propagator of eq. |Il] in 
the Schwinger Dyson equation for the quark self energy 



--V 



--V 



jk — mo 



(14) 



w.k w,k & D 

Another equivalent method is to use the Hamiltonian for- 
malism for the quark and antiquark creators and annihi- 
lators jlO]], and find the Bogoliubov Valatin transforma- 
tion which would minimize the vacuum energy density. 



In that Hamiltonian formalism the mass gap equation 
is also obtained when the quark antiquark pair creation 
operators are postulated to vanish in the Hamiltonian, 
in order to ensure the vacuum stability against sponta- 
neous generation of scalars. With the present method we 
project eq(|l4|) with the spinors ■ ■ ■ u and • • • v, and 
we get directly the quark and antiquark energy and the 
mass gap equation, 

E(k) = ul(k) \kk-a + m P - / - ,„ *yr iVi(k - k') 



2tt (2tt) 



^A-(fc')^/ 



^A+(fc / )r> ; 

w' - E(k') + it -w' - E(k') + it 

dw' d 3 k' 



= u\(k) |fcfc • a + mo/3 
fi|A+(fc')fli 



2tt (2tt) 3 
^A-(fc')^/ 



u s (k), 
iV(k-k') 
v s »(k), (15) 



_w'-E(k') + it -w' - E(k') + ie 

In the case of an instantaneous interaction, the loop inte- 
gral in the energy w removes the pole in the propagator 

dw i 1 

2tt w - E{k) + it ~ 2 

and in the case of a quadratic interaction, the loop inte- 
gral in the momentum is transformed in a Laplacian, see 
eq.ffl). Some useful properties are, 



(16) 



u\u s > 



1S SS > , ulv s ' 







a ■ ki<j2 



ulpv s , = -C 



, u\a ■ kv s > = s 



u\f3u s i = SS. 
u\d ■ ku s i = CS, 

u\f3d ■ ku s i = 0S ss i , u\(3d ■ kv s i = 1 
We get finally for the quark energy, 



a ■ kiu2 
a ■ ki<J2 
a ■ kia2 



(17) 



E(k) = kC + m S- 



1 



-SA(S) + Ck ■ A(kC) 



U 
~2 



1 r 



-35A(S) - Ck ■ A(kC) 



U ip 2 C 2 

-^ + kC + m S - — - — 



SCA^ ,, 



C 2 
k 2 



s 2 + lir 



(18) 



where in color singlets the U/2 term will cancel the — U 
term from the 2-body quark potential. For the mass gap 
equation we get, 



= \kS-m C + l 



-CA{S) + Sk ■ A(jfcC) 



1 



-3CA(S)-Sk-A(kC) } 
2SC 



a ■ kia2 



-Atp + 2kS- 2m C 



k 2 



2C Z + 1) A^ + 2SC ip 



(19) 
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The mass gap equation is in general a nonlinear integral 
equation, but in this case of a harmonic potential it sim- 
plifies to a differential equation. We solve it numerically 
with the Runge-Kutta and shooting method, see Fig. || 
for the solution. 




FIG. 5. We show the BCS chiral angle in units of Ko = 1. 
We also represent with a doted line the chiral angle that we 
obtain going beyond BCS. 



D. The pseudoscalar and scalar solutions to the 
Salpeter equation 

The homogeneous Salpeter equation for a meson (a 
color singlet quark antiquark boundstate) is, according 
to Appendix C, 

+ M(P)-E{k l )-E(k 2 ) ^ +{kp) = _i u t( fel ) x (jfe j p) 1 ,( fe2 ) 
-M(P)-E( kl )-E( k2 ) (f) -t^ p) = _ ro t (fcl)x(fcj p )u(fc2 ) 



X(k, P) = $& Vi(k - k')Q t {u{k[)cp + {k\ fV(*a) 



+v(k' 1 ) ( f > - f (k',p)u\k' 2 ))n l 



(20) 



where k\ — k + -j , k 2 — k — and P is the total mo- 
mentum of the meson. We use the Bethe-Salpeter am- 
plitude x as an intermediate step to compute the contri- 
bution of interaction V to the boundstate equation. The 
wave functions <f) + and 4>~ are equivalent to the Bcthc- 
Salpeter amplitude \ . For color singlets the contribu- 
tion of the infinite infrared constant U cancels, see Ap- 
pendix A. The equation is also flavor independent, and 
we will now concentrate on the momentum ® spin part 
of the wave- functions. We will now drop the U term and 
the color dependence from the equations. In this section 
the matrices f2; will only include the Dirac structure, 
f2; (g) tti = 70 ® 7o + 07 • (817. With the aim of studying the 
fo decay in a pair of 7r, we will now solve the bound state 
equation equation for the scalar fo in its center of mass 
frame and the equation for the pseudoscalar groundstate 
7r in the limit of small P and in the limit of large P. 

Due to the large mass of the scalar meson fo in this 
model, it turns out that the negative energy $r compo- 
nent for the groundstate is less than 10% of the positive 



energy component <j) + . The Schrodinger limit, where only 
the positive energy component is considered, is therefore 
acceptable. A general form for the 3 Po wave-function for 
the scalar is 



(f> + {k) SlS2 = k(j> s (k) 



[a ■ fctggjait 

V2 



(21) 



the truncated BS amplitude is then, 



2k(j} s (k) 
1 + S 



V2 



\/l 



we get for left hand side of eq. ( |2fj| ) 

u\k) x (k,0)v(k) = ^AU- 



1 



a 



1 - 3a 



CA (C0+) 



2 

1 -3a 



(22) 

5 A (S<p+) 
C 2 4> + (23) 



and the radial Salpeter equation for the scalar in the 
center of mass is, 



2E(k) - M 



dk 2 



k 2 ~~T + ~k 2 ~ 

1 + 2C 2 ., 



dk 2 



k 2 (p s = (24) 



Solving the bound state equation we find that the solu- 
tion of the equation is very close to a Gaussian, 



, K- 1 = 



5/2 ' 



.476 



(25) 



and the mass is M = 2.94 Kq — 970 MeV which is close 
to the most probable experimental mass of the /o ground- 
state. 

We now study the pseudoscalar groundstate in the low 
P limit, which was already studied extensively in the 
literature [|ll],12|. For vanishing P we find that 4> + — 
—4>~~ and both are proportional to sin(tp). This is due 
to the Goldstone boson nature of the n, see the result 
of Appendix C. However this component of the wave- 
function has zero norm, and it is necessary to include the 
next order of the expansion in P to determine the norm. 
The most general low P pseudoscalar wave-function is 
then, 



= Af- 1 



S- 



<r=K 



, M ( p ) t ■ P t S\i*2 

fi - m— • fc x a 1 



k J1 "^k / V2' 

where the norm is a function of the 7r energy, 



(26) 
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W 2 = (2f^) 2 M , M 2 (P) = M 2 (0) + P 2 



M 2 (0) = - 



2m (ipil)) 



ft 



d 3 k 

(2tt)3 



5 (27) 



and where in the case of an instantaneous interaction 
there are P,@ usually 2 different and fi s \ 



J 71 



(t) - 



(28) 



,(s) ,(t) _ i - 

^ ~ V 27T 3 



dfc -jfe 2 S<p + 4M7(flr 1 + l/2). 



Substituting the wave- functions of eq. (|26|) in eq. (20), 
and expanding the resulting equation up to the first order 
in P, we get the equations for the f\ and g\ components, 



h 



-kS+(2Ck)fi 

~i2 j i.2 i o2 * 2 i -2 



5i = 



l + a(2S' 2 + l) 

+4a (-SCAip - C 2 /k 2 + S*V + (p 2 j2) h 

-AaSC<p(f\ - h/k)] , (29) 
1 

l + a(2S' 2 -l) 
+2a{S 2 /k 2 - 2SCAtf - 2S 2 ip 2 ) 9l 

-AaSOpig! - gi/k)] . (30) 



kC+ (2kC + 2S 2 /k 2 ) gi 



It turns out that the parameter a has little effect on most 
functions, except for f\. The homogeneous equation for 
fi has the solution ao = —.195 an thus fial/a — do- This 
will essentially affect fi , f„ and the pion velocity c. 



We find for a ~ -.18 that c = 1, fi t] = f, 



0.21^ 



69 MeV. This shows a clear improvement of the model, 
with a correct relativistic pion and a better f v . 




FIG. 6. We represent the n wave-functions functions S, 
fi/k and gi/k respectively with solid, dashed and dotted 
lines, in the adimensional units Ko = 1. 



We now discuss the pseudoscalar groundstate in the 
other limit of large momentum P. In this case the neg- 
ative energy components are suppressed by a factor of 
1/P. The chiral angle <p, depicted in Fig. || vanishes 
completely, and the spinors are simpler, for instance, 



u s (ki) 



1 + a ■ ki 
V2 



it s (0) , ki 



P+pki 



(31) 



where the index ±_ denotes the projection k — (k • P)P of 
a vector k in the plane perpendicular to P. The vertices, 
up to first order in 1/P are for instance, 



1 



ul(ki)u s > (k[) ~ <5 SS ' - —ia ss >P x (kx - k'j_) 



P 



(32) 



4(fc 1 )au s /(fc 1 )~P^ 



(fcj_ + k'±) + ia x (fcj 
P 



k' 



Up to first order in 2/P the equation for positive and 
negative energy boundstate functions <f> (k)i<T2 / i s > 



2fc^ 
P 



T M-(l-a)A k )<j> ± (k) 



P 



-2i(l-o)p x Vfc-ja,^^)} 



-aA k (f> T (k) - -p^ffx ■ (j) T (k)a± 



(33) 



We find that the wave-function has a component with 
structure ia ■ P x k. However this component is smaller 
than the s-wave component by a factor of less than 1/P. 
Thus we find that the momentum (g> spin solution up to 
highest order in, is essentially a positive energy Gaussian 



function 4> p {k) 



v'2 



-P (34) 



For large momentum P we find that cf> + is quite flat in 
k, while 4>~ is almost negligible, 



2p 



A^ 



(35) 



This result is consistent with the relativistic space con- 
traction. We checked that the components that we ne- 
glected here would yield a small contribution to the /o 
decay. 



E. The coupling of a scalar to a pair of pseudoscalars 

The form factor P(P) for the coupling of a scalar /o to 
a pair of 7r can be decomposed in diagrams where a quark 
(antiquark) line either emits (absorbs) a pseudoscalar or 
a scalar. We use the truncated Bethe-Salpeter amplitude 
X, as an intermediate step to compute the coupling of a 
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meson to a quark line u^x u - F(P) is represented with 
a large triangle and the boundstate amplitudes \ an d 4> 
are represented with small triangles, 




(36) 

This includes the qq pair creation or annihilation of Fig. 
^. The same irreducible interaction for quarks which is 
used in the boundstate equations is also J2(| used for 
the annihilation. In eq. (pq ) the loop energies are triv- 
ial, see for instance eq. ( |A3| ), and we now compute the 
momentum ® spin contribution. 

We first consider the limit of a coupling to pions of low 
momentum P. The coupling of the pseudoscalar to the 
quark is derivative and thus is suppressed. For instance 
in the case of a massless 7r in the center of mass, using 
the wave- function of eq. ((261), we find that 



X 



Sa^Af- 1 A (5)1 /3 7 5 + o(P) 



u\k) X (k,0)u(k)=o(P) 



(37) 



which is consistent with the derivative coupling of a pion 
to a quark. The dominant contribution includes the cou- 
pling of the 3 Pq scalar meson to the quark (antiquark) 
line. The coupling F(P) to a pair of pions with low mo- 
mentum P is, 



F = tr 



d 3 k 
(2^)3 



4 (u f XsUq 



(38) 



where for instance we get for the scalar coupling 
tv (k)xs(k)u(k) to the quark line, 



2^2 



(1 - a)Ck ■ A(kSk<j) s ) 
-(l-3a)SA(Ck<f> s )] (39) 



except for the —i factor which goes with any potential 
insertion according to the Feynman rules. We will discard 
it in this section. Integrating by parts the eq. (|38| ) can be 
simplified with the help of the mass gap equation and we 
get, 



F = 



d 3 k k(f) s 
(2^72 



-kS 



dk 



aSCA 



(VV) 



O.Ollifo 2 - 0.052A/ 2 (P) 1/2 
/|M(P) 



(40) 



This coupling is very sensitive to the pion decay constant 
/„■ and to the energy M{P) of the pion. 



We now consider the opposite limit of large pion mo- 
mentum P. In this case the negative energy component 
of the pion is quite small. We expect that the dominant 
diagrams are the ones of the first line of eq. ([mJ), which 
include only the positive energy (j> + . This is present in 
either the coupling of a pion to the quark line, or the 
coupling of a 7r to an antiquark line, 



F(P) = 



d 3 k 
(2^)3 



Mh, o)(-A' fc ) 4> p {k', P)Mk, -p)U 



+^ p (k,p)<p p (k',-p)tq 



(41) 



k'=k 



where i g and tq are respectively the traces in Dirac in- 
dices, 



tq — trl — j=-v ' kiu^ (ki)Qiu(k' 1 )— -^v\k' 2 )Sl.iu(k2)— -i= 



i a 2 



i cr 2 



V2 



V2 



= tr< — = = ih = 

\ V2 2 V2 V2 

1 (1+/3) {l + a-h) 1 



V2 2 



-75- 



V2 



V2 V2 



V2 



(42) 



and tq, of the second diagram , with the coupling of the 
pion to the antiquark yields the same result. The total 
coupling is a functional of the scalar and pseudoscalar 
wavefunctions, which are described in eqs. by 
Gaussians. We now apply the Laplacian to the functions 
of k' . The dominant term of the expansion in 2/P comes 
from the derivatives of 4>p. A derivative of (f>p(k',P) 
will be proportional to P/2 when the Gaussian integral 
is performed , while a derivative of t q or tq which are 
functions of respectively k\ or k' 2 is proportional to 1/P 
and will not produce a dominant term. The traces then 
simplify to, 



U 



V2 



(43) 



We now apply the Laplacian to (j)p, and expand the k\ ■ k 2 
in a series of 1/P. It is convenient to define = 2a 2 . + 
a 2 , and we get finally find for the momentum £g> spin 
contribution, 



F(P) = 



-647T4 aial P 



„3 2 



3^ « 



(44) 



where the dominant term is the a term which is of the 
order of 0.43 K~ 1/2 for P of the order of 2K . 

The color factor for/ and tt color singlets is 
The flavor factor for the coupling of a scalar isosin- 
glet (uu +_dd)/y/2 to a pair of pseudoscalar isovec- 
tors, say ud and —du, with a flavor independent quark- 
antiquark annihilation is — 1/ The total coupling is 
then -F(P)/\/6. 
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The function F(P) is very cumbersome to derive in the 
case of intermediate momenta. For momenta of the order 
of Kq, see Fig. f7[ matching the high and low P limits 
with an interpolating function is a possible approxima- 




FIG. 7. We show F(P), in the adimensional units of 
Kq = 1. The dotted line and dashed line correspond respec- 
tively low P and high P limits. 



F. The / (980) -> tttt decay 

The decay width of a fo in a pair of 7r can be calculated 
from the the Breit-Wigner pole in the meson propagator. 
We c all bare the ladder meson, see eqs. (0), ( pl[ ) and 
( |C13| ). The bare mass Mo is real and is a solution of 
the Bethe-Salpeter equation. When coupled channels of 
mesons are included, the bare meson is dressed. The 
dressed pole is composed by the bare mass Mo plus the 
coupled channel contribution which includes a real mass 
shift, and an imaginary term in the case where the mass 
is above the coupled thresholds. The mass is then, 



M = M + AM , AM 




where we only included a loop of bare 7r, which is the sim- 
plest contribution to the scalar self energy. The integral 
in the loop energy provides an extra 2iri factor, which 
verifies that in general AM includes a real component, 
and we find that. 



AM = 6 



EM! £M 



(27r) 3 M /o -2^(7 2 + M2+i € 



(46) 



where the factor of 6 includes the 3 different flavors of 
the isovector 7r, and the factor 2 from the direct and 
exchange diagram of the self energy in eq. (|4^) . In this 
section we represent as usually the width of a resonance 
by a r, which should not to be confused with the same 



symbol which is used in other sections for the vector or 
axial vertices. The width T is a simple function of the 
imaginary component of AM, 



Im(AM) = -i- , T 



±-PM fa \F{Pf 



(47) 



where P = \ 



AM% is the momentum of the emit- 
ted 7r. Let us consider the case of a scalar mass Mf of 
the order of 1 GeV, where P is larger than the scale of 
the interaction by a factor of 1.4. In this case, it is sensi- 
ble to use the limit of large momentum P for P(P), see 
eq. (|44|). We finally find a partial decay width of just 
40 MeV for Fq — > nn which lies within the experimental 
limits. We expect that a complete calculation without 
the large P approximation would not deviate from this 
by more than a factor of 2. 

This is also compatible with the narrow resonances 
/o(1500) and ao(980) which are possible groundstates. 
Concerning Kq (1430) which is wider and has decay prod- 
ucts with a larger momentum P, the function r(P) of 
eq. (|47]) has the correct qualitative behavior of being pro- 
portional to P 3 for intermediate momenta. However in 
this model the exponential decrease is too strong, and the 
model needs some improvement in order to reproduce the 
correct #£(1430). 



III. GOING BEYOND BCS WITH FINITE 
COUPLED CHANNEL EFFECTS 

A. The Mass Gap Equation and the self energy 

We find in Appendix D that the minimal extension of 
the mass gap equation beyond BCS is achieved with a 
new tadpole term in the self energy, 



0(k,w) 




k — k 



(48) 



where the sub-diagram O is defined as an intermediate 
step. This amounts to extend the MCE for the self energy 
of the quarks with the simple one meson exchange. Using 
the Weyl fermions, and expanding the ladder in meson 
poles, we find that the self energy of the quark (anti- 
quark) has a diagonal component which contributes 
to the dynamical mass of the quark (antiquark), 



(49) 







and the energy of 1 quark is identical to the BCS one 
of eq. except for the expected changes of the chiral 
angle ip. In eq.(]49|) we only included the nonvanishing 
diagrams which remain from an expansion in powers of 
1/U. The free Green functions are proportional to U^ 1 . 
The interactions without a pair creation or annihilation 
are proportional to the infinite infrared constant U, while 
the remaining interactions are finite. It turns out that 
the new coupled channel diagrams vanish. This happens 
because in the limit of an infinite —U the box diagrams 
in eq. ( |48| ) which contribute to the quark energy vanish. 
This is the case for instance of diagrams (a) (c), 



0, 



0, 



(a) 



(b) 



(50) 



Eq — Erf remains the 



We find that the quark energy E 
BCS one of eqs. @ to @, 

i f r rfiy 
E{k) = u\k)a ■ ku(k) - -[^(fc) J j^V(k - k') 



A+(fc') - A-(fc') fiiit(fc)} 



(51) 



However the last diagram of eq. (|50|), which contributes 
to the mass gap equation is finite. 

The mass gap equation is obtained when we impose 
that the antidiagonal components S a of the self energy 
must cancel. As in eq.(15) this component in obtained 
with the projection of the spinors w and v. This pro- 
duces a function E a with the quantum numbers of a 
scalar, see eqs. ( |T9|) and (|2l|). In order to use the re- 
sults of the preceding section, it is convenient to fold E Q 
with a generic scalar wave function . Then the result- 
ing product must vanish for any <^ . In fact this ensures 
vacuum stability since this prevents the vacuum to decay 
in scalar modes. The diagrams that contribute to the an- 
tidiagonal component of the self energy in the mass gap 
equation are now, 



tr{^+(fc) f S a } 





(52) 



In eq.(|52j) we only show the diagrams which are nonva- 
nishing in orders of 1/U, and in fact they all are finite, 



of order U°. The first pair of diagrams are BCS dia- 
grams. It turns out that the new diagrams are the same 
diagrams which contribute to the /o — > 7r ® ir coupling, 
except for the negative energy wave function of the 7r 
and for the integral in the 7r momentum P. The neg- 
ative energy wave-function cj>~ always vanishes for high 
momentum P and in the case of low momentum P it 
is relevant only and for the pseudoscalar family of the tt. 
We suppose that the large number of excited states is not 
sufficient to compensate the smallness of <j)~ , and we will 
not consider this ultraviolet problem. Thus we will only 
include the coupled channel contribution of the 7r family. 
In the case of low momentum, the 7r family has a ex- 
tremely large cf>~ , however the coupling to a quark u*x u 
is derivative and vanishes. This prompts us to neglect 
the last line of eq. (|52"|) . Finally we can remove, with a 
functional derivative, the variational scalar wave function 
(f)~f o . The result after integrating in all the loop energies 
can be represented, 




where the lines only represent the spinors u, v, or 

v \ the integrals and the traces and no longer include the 
quark or anti-quark propagators. The mass gap equation 
= S ~ 1 a - S a is now, 

= +ul 1 (k)a ■ kv s& (k) 



+u 



,.(k')(- 



d 3 p 



2 J (2tt) 3 ^ S2S3 



(P,k'-P/2) 



,(fc')(% 



d 3 P 



(2tt) 3 ^ S2S3 

0-. 4 (P,*'-P/2) (54) 
where the sum over repeated spin indexes Sj is assumed. 



B. Model independent effects of the coupled 
channels 

The dominant effect of coupled channels is to multiply 
the potential term in the mass gap equation by a factor 



10 



of . 



d 3 p 

(2^)3^ 



-t 



(55) 



This clearly decreases the term which is the source for 
the spontaneous breaking of chiral symmetry. Thus the 
coupled channel effect is to restore partially the chiral 
symmetry. This effect is independent of the quark-quark 
interaction. 

The signs of the new terms in the mass gap equation 
deserve a special attention since they determine whether 
the coupled channel effect will increase or decrease the 
chiral condensation. Because the coupled channel terms 
introduce in the mass gap equation a new fermion loop 
it is natural for Dirac Fermions that the coupled channel 
terms should be affected with a — sign. 

When the Dirac fermions are translated into Weyl 
fermions the quarks divide into the species of quarks 
and antiquarks which have independent field operators 
and propagators, and the — signs are transfered from 
the propagator and the loops into the antiquark vertex 
and the exchange diagrams, and loops with quarks (an- 
tiquarks only). In this case we check with Weyl fermions 
that the — sign persists and is due to the quark (anti- 
quark) exchange. Only retardation, which was not in- 
cluded here, might perhaps oppose to this negative sign. 

This sign can also be understood from the perspec- 
tive of the Mexican hat potential — Act 2 + /icr 4 of effective 
meson models. In this case the quadratic term spon- 
taneously creates a scalar condensate, while the quartic 
term opposes to the condensation and the actual conden- 
sate corresponds to the minimum of the energy density 
where the two terms are balanced. In the present paper 
there are three terms, a kinetic term which opposes to the 
condensate (and has no correspondence in the effective 
meson models) a BCS term which spontaneously breaks 
chiral symmetry (it is equivalent to the quadratic term 
of effective meson potentials) and a beyond BCS term 
which is equivalent to the quartic term in effective me- 
son potentials. This correspondence, which is supported 
by the mean field theory where (ipif^) — <J 1 confirms the 
negative sign of the coupled channel term. Thus we may 
assume quite generally that coupled channels oppose to 
the breaking of chiral symmetry. 

An interesting feedback from chiral symmetry to the 
narrow width the groundstate occurs. Chiral symmetry 
breaking can be understood variationally, the solution 
ip(k) of the mass gap equation also minimizes the energy 
density £ of the vacuum. £ is the sum of three terms, the 
free one, the BCS one and the beyond BCS one. Only 
the BCS term is negative and drives f(k) away from the 
trivial vanishing solution. The actual solution ip(k) min- 
imizes the free term and minimizes the beyond BCS term 
and at the same token produces the most negative BCS 
term. In eq.(|5^), we saw that the beyond BCS diagrams 
in the mass gap equation are similar to the diagrams of 
coupling of a scalar to a pair of pseudoscalars (with low 



to moderate momentum). Thus we conclude that this 
coupling is naturally suppressed and that this suppres- 
sion is selective in the sense that it is not supposed to 
occur in other hadronic couplings. This has an amplified 
effect in the scalar width which is a function of the square 
of this coupling. 



C. Solution of the mass gap equation 

We will now focus on the dominant terms among 
the coupled channel contributions, We obtain the 
momentum (g> spin coupled channel contribution, 

f d Pd W - 2P2 ^ k + P JT ik + m ' (56) 



(27T) 2 



From the eqs. (26 34|) we get for the the integrand of eq. 



P 

P - 



0.14 



S(ki)-M(P)£j^ 

a 2 p a/4 fc?/ ; 



_ i q - - P 

L - J (l-a) 11 /4 e 



(57) 



where the spin factor ^| —ntf — \ is included. The inte- 
gral that leads to £ is now evaluated with an interpolation 
between the two limits of eq(^7|), see Fig. |[ Because this 
interpolation is arbitrary, we have to include in the result 
a theoretical error. 




FIG. 8. We show the integrand of f(0), in the adimen- 
sional units of Ko = 1. The dotted line and dashed line 
correspond respectively to the cases where cf> is obtained in 
the low P limit and in the high P limit. 

Wc obtain the momentum® spin coupled channel con- 
tribution, 



f(fc) = -(0.4 ±0.1)5 



(58) 



which turns out to have a shape very close to the func- 
tion S = sin((f>) which was evaluated at the BCS level. 
We estimate that the coupling to the quasi Goldstonc 
bosons, including the momentum, spin, color and flavor 
contributions yields, 



ff 



£(fc)^ -(0.4 ±0.1) S± (Nj-jjr 
~-(0.3±0.2)S, iV/ = 2^3 



(59) 



Nf is the number of almost massless quark flavors which 
empirically is between 2 and 3. 

The mass gap equation can be solved for a coupled 
channels contribution equal to £, when £ > — 1, 



= u f 



a-k-Uli j I '(A" - A id <)<>/ 



(60) 




(1+0+^ 



(l + O + aSCAcjo'-jfcio-a 

The solution y> for £ = —0.3 Sbcs is shown in Fig. ||. In 
this case and for the same parameter K$, quark conden- 
sate ('ipip) is decreased by a factor of (1 — 0.12) 3 . 



IV. RESULTS AND CONCLUSION 



the mass gap equation from coupled channels is closely re- 
lated to the coupling of a scalar to a pair of pseudoscalars. 
The feedback enforces that the width of a groundstate 
scalar decaying to a pseudoscalar pair (with low or mod- 
erate momentum) is reduced when compared to the width 
of any other resonance. This effect is model independent 
and contributes to understand the scalar meson puzzle. 

Concerning real mass shifts we estimate that they are 
canceled due to the new terms which are introduced by 
the Ward Identities. This might improve previous [fl||l2|] 
coupled channel calculations where this cancellation was 
not explicitly included. An systematic shift of the hadron 
spectrum is not expected. In this sense we agree with the 
results of Geiger and Isgur Q . Nevertheless mass split- 
tings between states with different quantum numbers are 
still expected from the boundstate equation. Different 
theoretical problems that could be reviewed with these 
new techniques are the contribution of the coupled chan- 
nels to the hadron spectra, for instance to the interesting 
Tj' mass or to the N mass. 

We conclude that in general the results of this paper, 
without ruling out other perspectives, explain why the 
simpler quark model is so successful. 



We developed a general formalism to include both the 
effects of chiral symmetry breaking and strong hadron- 
hadron interactions in quark models. This is encouraging 
since both effects are firmly established in phenomenol- 
ogy. We find new general effects in the scalar meson 
width, in the breaking of chiral symmetry and in the 
mass shifts of the hadron spectrum. Quantitative results 
are computed within a model which belongs to a class of 
Nambu and Jona Lasinio absolutely confining instanta- 
neous interactions, in the case where the coupled pair of 
mesons are accounted as bare mesons. 

We find in this model that the mass and width of 
the light qq scalar /o meson are close O] to the exper- 
imental mass and width of the /o(980), and not to the 
/ (400-1200) or the / (1370). This apparently indicates 
a possible solution to the scalar meson puzzle without 
meson molecules, glueballs or strongly nonlinear coupled 
channel effects. In this case the attraction which is visible 
in 7r 7r phase shifts and in the intermediate range NN in- 
teractions would need other interpretations [^7j than the 
very wide a meson. 

Compared with x^B at the BCS level, a new param- 
eter has been identified, which leads to the percentage 
of coupled channel effects in the mass gap equation. We 
find that coupled channels suppresses the breaking of chi- 
ral symmetry. This results are model independent. With 
our model we get a suppression of the quark condensate 
by 5% — > 55% when the coupled channel effects are in- 
cluded. 

We find a new interesting feedback mechanism from 
chiral symmetry to coupled channels and explain it varia- 
tionally. The chiral symmetry restorating contribution to 
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APPENDIX A: CONFINEMENT WITH 
INFRARED FINITE COUPLED CHANNELS 

At the BCS level, the boundstates are obtained with 
the ladder approximation which is equivalent to the 
Salpeter equation, and to the Schrodinger equation (see 
Appendix C). In this case the infinite infrared divergent 
constant U of the interaction ([||]) is extremely conve- 
nient |^2| to remove the colored states, which have masses 
proportional to U. Let us consider the dominant terms in 
orders of U of the energy of a system of n quarks and an- 
tiquarks. The 1-body energy includes the self energy (|is| ) 
which is calculated with the Schwinger Dyson equation, 

-2-uIL-K, (Ai) 
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where A is the Gell Mann matrix of the total color of 



the system. The energy ( Al) vanishes for color singlets 
only, thus colored states have an infinite energy and are 
confined. 

To include coupled channels in the energy of a color 
singlet, for instance in a meson, we consider the complete 
series of diagrams that contribute to the irreducible q q 
interaction. One has to include all the possible number of 
quark loops and all the possible insertions of the micro- 
scopic quark-quark interaction. Then this series can be 
resumed in order to factorize the bare meson (q q) and 
hadron (qqq) ladders. According to Appendix C, the 
ladder needs integrals / ^ in all the external relative 
energies in order to have a hadron pole. To ensure in a 
particular diagram that the meson poles are present it is 
convenient to decompose the ladder in, 



+ 



(A2) 



where in the right hand side the ladder is integrated and 
contains the meson pole, and the first pair of diagrams 
contributes to the overlap interactions of hadrons. In or- 
der to simplify the calculations it is convenient to trun- 
cate eventually the series of diagrams. The pertubative 
parameter is then the number of considered ladders. This 
is both straightforward and also close to the hadronic 
phenomenology. 

We now show how the interactions of color singlet lad- 
ders can be finite, when they are built from the infrared 
divergent quark-quark microscopic interaction (E), and 
the quark propagators are divergent as well. When a 
hadron is emitted or absorbed we have quark-antiquark 
annihilation overlaps, for instance in a 3-meson vertex, 




(A3) 

We first integrate the relative energies in the first dia- 
grams of q. JA5|) , 



dwdw' 



(2tt) 2 w - E x + ie 14,3 -w' - E' 4 +ie 



w' — E' 3 + ie w — £3 + —w — £2 + ie 



- i Vi 4 , 3 G a Go b G c 



(A4) 



where the Go n are in fact part of the respective BS am- 
plitudes which are untruncated. The remaining factor 
—iVi4,3 is finite because the term proportional to U in 
the quark-antiquark annihilation vertex is vl(k)5 3 (k — 
k')u s >{k')=Q. 

When the number of hadrons are conserved, this is 
the case in elastic scattering, we have quark exchange 
overlaps, for instance in a 4- meson vertex, 





(A5) 



The integrals in the relative energies are, 
dw i i i 



2tt w — Ei + ie —w — £4 + ie w — £3 + ie —w — E2 + ie 
i(Ei + £ 2 + £3 + E^GoaGobGocGod 



dwdw' 



(2tt) 2 w-E 1 + ie 
i i 



' r H^3)-7 ' 



w' - E' x + ie 
i 



—w' — E' A + ie w' — E' 3 + ie w — £3 + ie —w — £2 + ie 
= —iVi3Go a Go b Go c Go d . (A6) 

The four Go„, will be absorbed by the Salpeter ampli- 
tudes and the remaining factors are, 

i{E x + £ 2 + £ 3 +£ 4 - ^13 - ^ 24 )Pi3 

= i{Ei + £ 2 + £3 + £4 + V u + V 23 )Pi 3 
= i£i 3 (£i + £2 + Vu + £3 + £4 + V34) 
= iP 13 (0E/ + ...) (A7) 

whe re th e infinite —U terms cancel in the same way as in 
eq. (Al) because the mesons 32 and 14 in the left and 12 
and 34 in the right are color singlets. We find that all the 
terms proportional to U cancel when the complete set 
of diagrams which contribute to the interaction between 
color singlets is included. 

In this framework the masses of bare hadrons and 
their interactions are finite and can be evaluated. Then 
we compute the masses and widths of dressed hadrons. 
These are the final freedom degrees which can be com- 
pared with the experimental spectrum of hadronic reso- 



APPENDIX B: THE LIGHT tt AND WARD 
IDENTITIES 

The solution to the pion mass problem is found using 
p0| , p9t the Ward identities (WI) in order to insure that 
the bound state equation for the pion -a Bethe Salpeter 
equation with coupled channels - is consistent with the 
non linear mass gap equation. The WI were first de- 
rived for fermion-gauge field theories, and were initially 

based on the simple observation that for free fermions 

1 

with propagator So(p) = — and a free vector ver- 



tex r£ = 7^ 



i{p„ - p')S(p)r»(p, P ')s(p') = s( P ) - s( P ') 



(Bl) 



The difference in the right member of the equation ex- 
tends the identity to renormalized propagators and ver- 
tices. This identity is then crucial for the conservation 



13 



of electric charge. The WI enforce that the self energy 
of the MGE is obtained (without double-counting) from 
the BS kernel by closing the fermion line where the ver- 
tex is inserted. Inversely, they also ensure that the BS 
kernel is obtained if one inserts the vertex in all possible 
propagators of the self energy. For instance this map- 
ping is trivial at the BCS level where the mass gap equa- 
tion (^) is clearly equivalent to the bound state equation 
(^). Let us now consider a more general case, where the 
fermion self energy include a product of bare propagators 
S ai i3 i (ki + P). The external momentum is P and fcj is a 
loop momentum. The propagators can be factorized and 
we get, 



S(P) = ■■■\{S ai p i (h+P) . 



(B2) 



Then the vertex T can be constructed if we insert a bare 
vertex in all possible bare propagators, and we get, 

r^Pi - p 2 ) = II S^Aitei + p i) 

3 il<J 

7 M II S aaf3i2 (k i2 + P 2 ) . (B3) 

i2>j 

Inversely the original self energy can be recovered if we 
substitute the bare vertex by the difference of propaga- 
tors, 

i(Pi f ,-P 2 ,)^=---J2 II S ail0il (ka + Pi) 

3 il<3 

[s aj0J + Pi)- s ajP] (k 3 + p 2 )) n s Q , 2te (k l2 + p 2 ) 

i2>j 

= [S(P 1 )-S(P 2 )] (B4) 

where the products which depend on both Pi and P 2 
cancel, and this removes the double counting. 

There is also a WI identity for the free axial vector 
vertex r^ J = 7^7 5 that involves the free pseudoscalar 
vertex = j 5 , 

-i{ Pll -^)S{ P )T^(p,p')S{ P ') 

+2imS(p)T 5 (p,p')S{p') = S(p)j 5 + -f 5 S{p') (B5) 

which is valid in a renormalization program providing 
the interaction is chiral inv ariant . In this case an equa- 
tion analogous with eq. ( B5 ) is complied by the axial and 
pseudoscalar vertices and by the self energy. 

A key product of the axial WI is the proof [Q,^o) that a 
pseudoscalar Goldstone boson exists when current quark 
are vanishing and chiral symmetry breaking occurs. The 
full propagator is then renormalized and the self-energy 
£ has a mass-like term, 

5- x (p) = S^ip) - E(p) , i£ = A(p)~ £>B{p) . (B6) 



S- 1 {p)=S Q \p)-V{p) = 



A(p) t>-B(p) 



(B7) 



If we substitute this propagator in the WI, we find the 
solution for the pseudoscalar vertex T 5 with a vanishing 

P-P 1 , 



r 5 ( P = p') = 



B(p) . 
to 



"7 



(B8) 



which diverges for a vanishing quark mass too and shows 
that the pole of a massless pseudoscalar meson appears 
in the axial vector vertex, with a boundstate truncated 
amplitude of 

P75 



X{p - p' ~ 0) 



(B9) 



where includes a norm. Incidentally, this identity also 
offers a proof of the Gell-Mann, Oakes and Renner rela- 
tion. If we expand the vertex T 5 in the neighborhood of 
the 7r pole, 



r 



X 



(p - p 1 ) 2 - M 2 



tr{ X S(p)r 5 S(p')} 



(BIO) 



where we included the integral in the trace, substituting 
T§ in cq. (B8) and performing a trace with S"f^S, we 
find for small masses 



M 75 



PS75S 



i f BS l5 S 
tr{ — ; 75 



= tr < 75 



BSj 5 S 



J —M 2 { J L TO 

=> -m tr{S} = -2m (ip V> = M 2 f 2 (Bll) 

In QCD it is necessary to include the axial anomaly, in 
the flavor singlet WI which corresponds to the 7/ channel. 
The flavored axial currents remain unchanged, in partic- 
ular the 7r, K, rj remain quasi Goldstone bosons. We will 
now simply assume that the 77' is heavier than 
the usual N 2 — 1 Goldstone bosons, where Nf is the num- 
ber of light quark flavors. 



APPENDIX C: THE SALPETER EQUATIONS IN 
THE ENERGY-SPIN FORMALISM 

When the interaction is instantaneous, a simplification 
occurs in the Bethe-Salpeter S matrix in the ladder ap- 
proximation, 




E 



-E 



(CI) 
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The S matrix only has 4 independent sub-matrices, which 
have to be calculated iteratively. The other 12 sub- 
matrices are directly computed from the independent 4 
ones. The Salpeter equations are obtained whe n all the 
relative quark-antiquark energies of the system (CI) are 
integrated, 



dw 



2tt w + W/2 -E q +iew+ W/2 -E q + ie 



W-E„ 



(C2) 



The Salpeter equations can then be written in the com- 
pact form, 



S = Go + G VS 
=> (Go 1 -V)S = 1 
(Wa 3 - H + ie)S = i 



(C3) 



where 03 is the Pauli matrix, 
3 



s- 



5 H 

s~ 



G 

V 
H 



W-E q -Eg+ie 





-W-E„-Eg+i<L 



JVd fV a 
JVa JV d 

E q + E q + J V d 

JVa 



E„ 



JVa 



Eq +JV d 



(C4) 



and it turns out that H is an hermitian and positive 
operator. 

The Salpeter wave-functions are the solutions of the 
homogeneous equations, 



(Ma 3 - H)4> = , . 



(C5) 



which is an eigenvalue equation, similar to the Schro- 
dinger equation, except for the 03 . This formalism is 
known as the energy-spin formalism, where <p + is called 
the positive energy wave-function and (\>~ is called the 
nega tive energy wave-function. The Salpeter equation 
(C5) is equivalent to the variational equation, 



= , 



(C6) 



which suggests that a normalizing condition of the wave 
functions might be, 



v 3 i0> = i^ + i 2 -in 2 = i 



(C7) 



Let us study the class of solutions <p u where this norm 
is possible, i. e. |</>+| > |<^>~|. Once the single quark 
energies E and the 2-quark diagonal Vd and antidiagonal 
V a potentials are defined, the solutions can be obtained 



numerically, either iteratively or variationally, and one 
finds a whole spectrum of solutions with energy M > 0. 
When M increases, we find that \<f)~\ is proportional to 
M , and in the limit of large mass we exactly recover 
Schrodinger equation. However another class of solutions 
<j>d is unavoidable, with a one to one correspondence with 



the 



and where 



Mo"3 



> 



1 



(j>+ , M d = -M u 
(^d|cr 3 |0d) = -1 



(C8) 



thus the spectrum is unbound. Exactly half of the solu- 
tions have a negative mass and a neg ative norm. When 
we include the infinitesimal ie in eq. (|C5|), we find that 
the larger component of (j) dominates and the eigenvalues 



M u -» M u - ie 

M d -» M d + ie = -{M u - ie) 



(C9) 



The operators H and (73 are both hermitian, thus the 
set of solutions \4>) constitute a basis of the Hilbert space, 
orthogonal in the sense that ((j)\a 3 \(j>') = 0, and the iden- 
tity element is, 

i = ZS = £^><0|-£^><0| (cio) 



Inserting this partition in the S matrix equation ( jCJj ), 
we find, 

i = {Wa 3 -H)S 
= £ y—j— — 7— r (W^(73 -H + *e)a 3 ^ jjT^azS 



W3\ 



0-3 



0-3 



0-3 1 



*<0k3| 



W- M, 



(Cll) 



and the states with negative energy and negative norm 
can be reinterpreted as boundstates with positive mass 
moving forward in time, where the variable W = Pq in 
the propagator turns out to be negative, 



s=Ei40 

u 

+Ol\4> u ) 



W - M u + ie^ 
i 

-W -M u + ie { 



01 



(C12) 



and it suffices to work with the <p u . Diagrammatically we 
get for instance, 
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dwdw 



+ + 





(C13) 

where w and w' are the external energies. The quark 
antiquark BS amplitude is represented by the triangle 
— ►— D> and the boson propagator of the meson is repre- 
sented by the double line . It is consistent to sub- 
stitute the d wave-functions and masses in terms of the 
components <p+ , 4>~ and the masses M u . In the paper 
this reinterpretation will be assumed and we will skip the 
subindcx u . The practical result is that we can use the 
amplitude <fi + when the bra(ket) precedes (succeeds) the 
meson propagator, and we use the amplitude <p~ for the 
opposite case. 



APPENDIX D: THE BETHE SALPETER 
EQUATION WITH COUPLED CHANNELS 

We now go beyond BCS including the mesonic coupled 
channels both in the MGE and in the bound state equa- 
tions. It is convenient to return to the Dirac formalism 
in order to reduce the number of diagrams when we ap- 
ply the WI. We will extend the BS for the qq boundstate 
of the quarks with the minimal meson loop of coupled 
channels. 



.= = r + u_f= 




D 



(Dl) 



where we follow the approximation of including only 3- 
legged effective meson vertices in the meson interaction. 
We will show that this is imposed by the box at the left 
of the beyond BCS diagram in eq. (Dl). Upon iteration 
the equation, 



r = T + VG r + VG VMT 



(D2) 



can be resumed. If we factorize the S matrix at the ladder 



level So = (1 - G Q V)- X G Q , see eq. qCSJ), the ladder will 
appear in the middle of the coupled channel terms, 



G r = s r + GoV s v Mr 
r ~ v s Q r + v So v my. 



(D3) 



When the ladders So , including the two ones of the me- 
son loop in M, are expanded in meson poles and wave 
functions according to eq. (C12), we recover the meson 
- meson pair coupling of eq (|36[). The resulting pole of 
the vertex T is the mass of the dressed meson Mo + AM. 
This is equivalent to the Resonating Group Method equa- 
tions |26| ] for coupled channels of one meson with a pair 
of non- interacting mesons. 



We now use the WI prescription of removing the ver- 
tex and closing the respective fermion line to arrive at the 
mass gap equation ([l8|). To recover the full vertex equa- 
tion we must insert the vertex in all possible propagators 
of eq. (pig]). We then arrive at a WI consistent Bethe 
Salpeter equation for the vertex and for the boundstate, 



•= = r + + 



+ + ; " : • — : — ^ 



+ 




+ 



(D4) 

where the diagrams are shown in separate lines according 
to their properties. Line 1 corresponds to the BS equa- 
tion at the BCS level. Without the other lines it would 
reproduce the ladder Bethe-Salpeter equation for the ver- 
tex. The lines 2, 3 and 5 were separated because they 
contain all the contain terms proportional to the infrared 
divergent U, but they cancel in each line and all lines are 
finite. The remaining lines 4 and 6 contain the terms that 
one would except in coupled channel equations where a 
pair of mesons is created and then annihilated (except for 
the first diagram of line 4 which vanishes). With them 
we calculate for instance the partial decay width of a res- 
onance into a channel of 2 mesons. The last lines 5 and 
6 are only relevant for flavor singlets because the quark 
pair in the incoming meson is annihilated, thus for flavor 
vectors they are null. 

The cancellation of the infrared divergences become 
clear in the Goldstone-Weyl formalism. Let us consider 
for instance the diagrams of line 2 in eq. (D4), 




(D5) 



1G 



These 3 diagrams are infrared dive rgen t but their sum is 
finite, in an analogous way to eq. flAq) . 

It is important to remark that in the previous calcula- 
tions in the literature, the extrapolation from the ladder 
level to the coupled channel level would only include the 
diagrams of lines 4 and 6. We now find the previous 
choice arbitrary since the diagrams of lines 2, 3 and 5 
were not considered. The role of these diagrams is to 
cancel any real mass shift of the it due to the usual cou- 
pled channel diagrams of lines 4 and 6, in order that the 
7r remains a Goldstone boson in the chiral limit. This 



is ensured by the WI, see eq. (B8). Concerning widths, 
since the new diagrams are real, the results of the previ- 
ous calculations in the literature are correct. However we 
find no systematic real mass shift in the meson spectrum 
due to coupled channels. This contradicts most of the 
real mass shifts of hundreds of MeV which are common 
in the literature. Only splittings between different lev- 
els, due to neighboring cuts, may be affected by coupled 
channels. 
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